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Abstract 

Approximative analytic solutions of the Dirac equation in the ge- 
ometry of Schwarzschild black holes are derived obtaining information 
about the discrete energy levels and the asymptotic behavior of the 
energy eigenspinors. 

Pacs: 04.62. +v 

The problem of the Dirac perturbations in the gravitational field of a 
Schwarzschild black hole is less studied since the free Dirac equation in this 
geometry is not analytically solvable at least in the systems of coordinates and 
tetrad gauge fixings used until now. However, despite of this difficulty, inter- 
esting results were obtained using suitable analytical and numerical methods 

[l]-[3]. 

In this letter we would like to consider this problem in a new conjecture 
in which the local frames of the central static charts are denned by the tetrad 
fields we have proposed some time ago [4, 5]. This tetrad gauge shows off the 
central symmetry as a global one such that the separation of the spherical 
variables can be done as in special relativity [6]. After the separation of the 
angular variables we remain with a pair of simple radial equations depending 
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only on the metric components. In this approach we have found the energy 
eigenspinors of a massive Dirac field on the de-Sitter [7] and anti-de Sitter 
spacetimes [4, 8]. Unfortunately, this can not be achieved in the case of the 
Schwarzschild background. 

For this reason, our purpose is to find reasonable minor approximations 
of the radial equations in which these could be analytically solvable. In this 
way we obtain interesting information concerning the discrete energy levels 
and the asymptotic behavior of the energy eigenspinors. The results are 
presented in natural units with h = c = 1. 

In general, a manifold with central symmetry has a central static chart 
with spherical coordinates (t, r, 6, 0), covering the space domain D = D r xS 2 , 
i.e. r G D r while 9 and <fi cover the sphere S 2 . For defining the local frames 
we use the tetrad fields defined in Refs. [4, 5] whose components depend on 
three arbitrary functions of r, denoted by u, v and w, which give the line 
element as 

dr 2 r 2 



ds 2 = w 2 



dt 2 -^—- —(dO 2 + sin 2 6d<p 2 ) 



(1) 



We have shown [4] that in these circumstances the separation of the spher- 
ical variables of the free Dirac equation can be done as in the case of the 
central problems in Minkowski flat spacetime. Consequently, the Dirac field 
of mass m can be written as a linear combination of particular solutions of 
given energy. Those of positive frequency and energy E, 



U E j, mj , Kj (x) = U E ,j, mj , Kj (t,r,9,(f>) (2) 
"" ;[/ + M<,,W) + rMt,,(W)]e-' H 



rw(r) 3 / 2 

are particle-like energy eigenspinors expressed in terms of radial wave func- 
tions f ± and usual four-component angular spinors These spinors are 
orthogonal to each other being completely determined by the angular quan- 
tum numbers, j and rrij, and the value of Kj = ±(j + 1/2) [6, 9]. Moreover, 
they are normalized to unity with respect to their own angular scalar prod- 
uct. We note that the antiparticle-like energy eigenspinors can be obtained 
directly using the charge conjugation as in the flat case [8]. 

Thus the problem of the angular motion is completely solved. We remain 
with the radial wave functions which satisfy two radial equations that 
can be written in compact form as the eigenvalue problem HT = ET of the 
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radial Hamiltonian [4] 



mw 



d 




H 



(3) 



u 



d 
~dr 




—mw 



in the space of two-component vectors, T = \f + ,f \ T , equipped with the 
radial scalar product [4] 



This selects the 'good' radial wave functions, i.e. square integrable functions 
or tempered distributions, which enter in the structure of the particle-like 
energy eigenspinors (2). 

Let us consider now a Dirac particle of mass m freely moving in the 
gravitational field of a black hole of mass M with the Schwarzschild line 
element 



defined on the radial domain D r = (r , oo) where r = 2MG. Hereby we 
identify the functions 



that give the radial Hamiltonian (3). The resulting radial problem can not be 
analytically solved as it stays. Therefore, if we desire to avoid the numerical 
methods, we are forced to look for suitable approximations of the radial 
equations. 

In order to find an efficient method of approximation it is convenient to 
introduce the new dimensionless variable 




(4) 









and the notations 



fx = r m , e = r E , 



(8) 
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which allow us to rewrite the radial problem as 



= 0. 



(9) 



Moreover, from Eq.(4) we find that the radial scalar product in the new 
variable takes the form 



(J^i, T 2 ) = 2r J™ dx [x + T\T 2 



(10) 



Near singularity, for small values of x ~ 0, the exact solutions can be ap- 
proximated by JF S = |/+, f~\ T which satisfy the approximative radial equa- 



tions, 



x 



1 d + K 



1 d + K 



/.-(*) 



o, 



(ii) 



where we neglected the terms of the order 0{x). Performing the transforma- 
tion T s — > J- s — UJ- S — |/ s + , f~\ T , with the help of the unitary matrix 



1 i 
i 1 



we obtain the simpler system of equations 



(12) 



(13) 



This can be analytically solved in terms of Bessel functions Ji and Yi as [10] 



fj~(x) = Kj - ijj, (ci y/x Ji + (2i\jx) + c 2 \/x Y i+ (2i\jx)^ , 
f~(x) = yj Kj + ijj, (ciy/x Ji_ (2i\jx) + c 2 \/xYi_ (2i\ 

where C\ and c 2 are arbitrary constants and 

1 



Z± = 2ie ± - , Aj 



(14) 
(15) 

(16) 



We observe that near singularity these wave functions are regular but their 
phases are not determined since oc x T2zt + 0(x) for x — > 0. 

The asymptotic behavior, for very large values of x, can be studied ap- 
proximating \J\ + x 2 ~ x and obtaining thus an asymptotic radial problem 
for T a = \ fa, fa\ T - The asymptotic radial equations can be put in the form 



ld_ 

2~ dx 



i !Si 

~ x 



H\x + 



X 



-n 



i A. 

2~dx 



x 



J a 



X 



= 0, (17) 



where we introduced the new parameter 8 G [0, 1) which could play the role 
of a fit parameter. Of course, if we rigorously consider the approximation of 
the order 0(l/x) then we must take S — \. 

Since Eqs. (17) can be solved for any value of e, we have to look separately 
for solutions corresponding either to a discrete energy spectrum when e < /i 
or to a continuous one in the domain [pi, oo). For finding these solutions, 
we diagonalize the term proportional to x of the operator (17) using the 
transformation matrix 



T 



1 

2^ 



- e Vt 1 + e 



where we denote v = ^/ fi' 2 — e' 2 . After the transformation T a 
\ fai fa\ T we obtain the new equations 



(18) 

f a = TT a = 




vx 



J a 



e/x(<5 



J a 



(19) 



that can be analytically solved. 

The discrete quantum modes with discrete energy levels e < /i are gov- 
erned by radial wave functions which must be square integrable. After a 
little calculation we find that Eqs. (19) allow the particular solutions 

f+( x ) = N + x 2s >e-» x2 1 F 1 (s ] +q + l,2 Sj + l,2vx 2 ), (20) 

f-(x) = N-x 2 ^e-^\F 1 (s J +q,2s, + l,2ux 2 ), (21) 

where the confluent hypergeometric functions iiq [10] depend on the follow- 
ing real parameters 

^{5 - 1) 



k 2 + v 2 + fi 2 (5 2 - 1) , q = v + 



(22) 
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These solutions satisfy Eqs. (19) only if the normalization constants iV ± 
obey 

NZ = v±I ~ g) (23 ) 
N+ Kju - ne{5 - 1) ' V } 

We obtain thus the transformed radial wave functions (20) and (21) that 
have a good asymptotic behavior and are regular in x = 0. 

These solutions may be square integrable with respect to the scalar prod- 
uct (10) only when we impose the quantization condition Sj + q = —n, n = 
1, 2, 3, .... This leads to the equation 

Jtf + v 2 + fi 2 (5 2 -l) + v + ^ ~ 1} = -n , (24) 

which is of the third order in v and can be solved at anytime under Mathe- 
matica or Maple, obtaining thus the asymptotic energy spectrum. However, 
the analytic formula of these energy levels is too complicated to be written 
here. For this reason, we restrict ourselves to give some reasonable esti- 
mations. First we observe that the energies levels, E n j, must satisfy the 
conditions 

2m ^ <^f <l<^ + 4- (25) 

Kj + fi 2 (l — d) 2 m 2 jj, 

The last inequality prevents the Dirac particles to get near to the singularity 
where the black hole could capture them. For /i < 1 and small values of v 
we have to use the approximative formula 



E n<i ~ m 



1 _ 



(26) 



Moreover, in the case of \i <^ 1 and 5 = \ we recover the Newtonian result 

/ 1 n 2 \ l-e 2 G 2 M 2 m 3 , , 

En ^ m { 1 -8 (n + j + l) 2 ) =m -— tf— (27) 

for a particle of mass m = ji/r and angular momentum L — j ± | moving 
on an ellipsoidal trajectory of eccentricity 
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These asymptotical results, including the exact solutions of Eq. (24), 
must be considered prudently trying to verify and improve them using nu- 
merical methods. We note that only the massive fermions can be retained 
in Schwarzschild gravitational fields on discrete levels while the massless 
fermions can not do this. 

The scattering modes corresponding to the continuous energy spectrum 
e e [fi, oo) are described by particular solutions that behave as tempered 
distributions. In this case we introduce the real parameter v = \J e 1 — ji 2 
instead of v. Following the same procedure like in the previous case we 
look for the solutions of the transformed equations (19) where we replace 
v = %v. Hereby we obtain the transformed radial wave functions in terms of 
Whittaker functions M P)S and W PyS as [10] 

f+(x) = Ct-M p+ , S] (2ivx 2 )+Ct-W M (2iux 2 ) (29) 

X Ob 

f-(x) = Cr-M p _ !S .(2ii>x 2 )+C^-W p _ !Sj (2ii>x 2 ), (30) 

Ob Ob 

where we denote p± = =f| — q while the constants satisfy 

Ci = Qjsj - g) Cj_ = o , 31 x 

Ci KjV + ie/ji(5 — 1) ' C2 Kju + iefi(5 — 1) 

These solutions are useful for analyzing the scattering of the Dirac particles 
on black holes. 

Hence we have all the elements we need to write down the approximative 
form of the energy eigenspinors in all of the above discussed cases. Starting 
with the approximative form of the transformed radial wave functions f^ 1 we 
have to calculate the functions f ± which should give the particle-like energy 
eigenspinors (2). The antiparticle-like energy eigenspinors can be calculated 
using the charge conjugation as [8] 

VE,j, mj , Kj = (UE,j, mj ,K 3 ) C = C(U E,j,m J ,K j ) T , C = 27 2 7° . (32) 

We hope that the approximative solutions presented here should help one 
to improve the analytical or numerical methods for studying the Dirac field 
in the Schwarzschild background. 
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